Magnetohydrodynamic properties of incompressible Meissner 

fluids 



A. Maeyens^, J. Tempere^'^ 
^TFVS, Universiteit Antwerpen, Groenenhorgerlaan 171, B2020 Antwerpen, Belgium, and 
'^Lyman Laboratory of Physics, Harvard University, Cambridge MA02138, USA. 

(Dated: June 18, 2007) 

Abstract 

We consider a superconducting material that exists in the hquid state, more precisely, in which 
the Meissner-Ochsenfeld effect persists in the liquid state. First, we investigate how the shape 
of such a hypothetical Meissner liquid will adapt to accomodate for an applied external field. In 
particular, we analyse the case of a droplet of Meissner fluid, and compute the elongation of the 
droplet and its quadrupole frequency as a function of the applied field. Next, the influence of 
an applied field on the flow of the liquid is studied for the case of a surface wave. We derive 
the dispersion relation for surface waves on an incompressible Meissner fluid. We discuss some 
candidate realizations of the Meissner fluids and for the case of a superconducting colloid discuss 
which regime of wave lengths would be most affected by the Meissner effect. 
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I. INTRODUCTION 



The superconducting materials readily available to man are all solids. Both for con- 
ventional superconductors such as mercury and for high-temperature superconductors, the 
critical temperature lies well below the melting temperature. Nevertheless a crystalline 
substrate is no necessary prerequisite for superconductivity. Even for the case of phonon- 
mediated superconductivity, it can be argued that liquid metals also exhibit phonons 
Indeed, for the case of dense liquid metallic hydrogen, a superconducting phase has been 
predicted 

In this contribution, we investigate how a liquid superconductor would respond to applied 
magnetic fields. More precisely, we consider a Meissner liquid, since we are interested in the 
magnetic response rather than the electric response. A candidate realization of a Meissner 
fluid might be a suspension of superconducting particles. In Ref. 0, 4|, such a suspension 
of micron-sized superconducting cuprate particles in liquid nitrogen was investigated, and 
when those particles were coated with a layer of ice, acting as a surfactant preventing the 
particles from coagulating, a Meissner liquid state was reported. 

The basic property of the Meissner fluid is its capability to expell magnetic flux B from 
its bulk by supporting a persistent surface electric current Jg. In type II superconductors 
above the first critical field, the flux is incompletely expelled. Here, we consider type I 
superconductors or type II superconductors below the first critical field. In Ref. jSj], it was 
shown that when equilibrium is disrupted, B and relax back to (new) equilibrium values 
with a characteristic time of r = 10^^^ s. The time scale for the hydrodynamic motion of a 
fluid surface (as determined by the fluid's density and surface tension) is much slower, and 
this allows to separate the dynamics of the currents and magnetic fields from that of the 
fluid surface. Hence, for the different problems investigated in the present study, we assume 
that at each time the magnetic field has reached its equilibrium value for the given surface 
deformation of the fluid 6|]. 

We start in Sec. II with the study of the shape deformation of a droplet of incompressible 
Meissner liquid placed in a uniform external magnetic field. From the energetics of the 
optimal shape we calculate the quadrupole mode oscillation frequency. In Sec. HI, we 
investigate the higher-mode frequencies by neglecting curvature and modeling surface waves 
on a Meissner liquid. We derive the modification of the wave dispersion by a magnetic field 
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parallel with the surface. 



II. MEISSNER DROPLET 

In the absence of gravity, a fixed volume of liquid will form a spherical droplet to minimize 
surface tension. When a magnetic field is applied on a droplet of Meissner fluid, the fluid will 
magnetize in order to expell the magnetic flux. The energy associated with this magnetization 
is expected to be smallest for a cigar-shaped droplet with the axis of cylindrical symmetry 
along the magnetic field. However, deforming the droplet into a prolate spheroid will also 
increase the surface tension energy. The optimal shape of the droplet can be found by 
balancing the magnetic and surface tension energies. 

As possible shapes we will consider prolate spheroids with the major axis parallel to the 
applied magnetic field. The variational parameters determining the droplet shape are a, 
the equatorial radius or semiminor axis, and c, the semimajor axis. We denote the applied 
magnetic field by H^. The droplet responds to the external field by magnetizing. The 
magnetization field M is zero outside and nonzero inside the droplet. Since then V ■ M is 
nonzero at the droplet surface, magnetic charges are induced on the surface. These charges 
give rise to a demagnetization field H^. 

Inside the droplet, the demagnetization field is straightforwardly related to the magne- 
tization through H^" = — n(a, c)M where n(a, c) is the demagnetizing factor for a prolate 
spheroid 0]: 



n(a,c) =1 
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where a = ^ (? — is the focus distance. The total magnetic flux is 

B = /i(H„ + H, + M), (2) 

with [I the vacuum permeability. Since inside the Meissner fluid B has to be zero, we find 
that the required magnetization satisfies 

M = ^H„, (3) 

1 — n(a, c) 



in agreement with the result of Ref. 

To find the demagnetization field HJ^"* outside the droplet, we solve V x HJ^"* = (since 
there are no induced currents in the volume outside the droplet). This implies that the 



outside demagnetization field can be written as the gradient of a scalar magnetic potential 



H 



out_ 



—Vif). The boundary conditions are that ^ at infinity and that the tangential 



component of is continuous along the boundary of the droplet. Equivalently, we can use 
that the normal component of the magnetic flux is continuous accross the boundary. The 
natural coordinate system to express these boudary conditions are the prolate spheroidal 
coordinates {rj,6,(j)}. The droplet boundary is then defined by fixing 77 = 1]^. The equation 
for the scalar magnetic potential becomes |9(] 

1 



Q;^(sinh^ 7] + sin^ i 

d'^ih ^ dib 



d'^ip dip 
—r + cothr/- — h 

■■ 0. 



dr] 



The solution of this equation can be found by separation of variables, yielding 

00 

i^iv,^) = ^^„<5„(coshr7)P„(cos6'), 



(4) 



(5) 



n=0 



where P„ and Qn are the Legendre functions of the first and second kind, and the An are 
integration constants. The boundary condition H^^ at infinity has already been used. 
The second boundary condition is 



H 



out 



[1 - n{a,c)] M ■ e^l 



where is the unit vector in the 77 direction. This corresponds to 



(6) 



e^- (cosh T]) Pn {cos 9) 
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The demagnetisation field outside the droplet is then 



-In L/^ 



\/ sinh^ ?7 + sin^ 9 



X |cos6' 
+ sin0 



coth ?7 + In ( coth — ) sinh 77 



1 + coshi] log(tanh 



eg 



(7) 



(8) 



(9) 
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Fig. [T] depicts the various contributions to the magnetic flux: the magnetization, the de- 
magnetization fleld, and the magnetic flux. The magnetic flux remains tangential to the 
surface, flowing around the boundary of the droplet. 
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FIG. 1: (color online, 2 column) Panel (a) shows the magnetization, which is only present inside 
the spheroidal droplet (shaded area), and which is equal to —nUa- The demagnetisation field H,^ 
is shown in panel (b). The total magnetic fiux B = /u(Ha + + M) is shown in panel (c). All 
arrows in panels (a) and (c) are scaled by the same factor, those in panel (b) are scaled by twice 
that amount for visibility. 



To calculate the energy difference Em between the unmagnetized and the magnetized 
droplet of a given, flxed shape with shape parameters a and c, we use the thermodynamic 
potential U = f/ + H ■ B with U the internal energy, such that dlf = B ■ dH and 



EM{a,c) 



Ha 



Bfr) ■ dUir) 



dr. 



(10) 



Using the results (T?),© for the flelds, the expression simplifles to 

1 f^Hl 



EM{a,c) 



1 — n{a, c) 2 

where V = Ana'^c/S is the volume of the droplet. This term would favour, for a given volume, 
a more elongated spheroid. It is counteracted by the surface tension energy Es = crS with 
a the surface tension and S the surface of the spheroidal droplet. This evaluates to 

-1 



Esia, c) 
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1 — a^/c' 
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(11) 
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To find the optimal surface, tlie total energy Es + Em need to be minimized as a function 
of {a, c} with the constraint of constant volume. If tq is the radius of a spherical droplet 
with volume V, we can introduce dimensionless parameters d = a/r^ and c = c/rg. The 
constraint of constant volume then allows to eliminate one of the variational parameters since 
V = Aira'^c/S = 47rro/3 leads to c = a~^. Writing the total energy as = {Es+EM)/i'2iTarQ) 
then leads to 



E(d) 



7^ + d^ + 



sin"^(Vl - a^) 



a6 



1 - n{a) ' ' aVl 
where F is a dimensionless parameter expressing the applied magnetic field 
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and 



n a 



1 



n6 
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a6 
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(12) 



(13) 



(14) 



Note that 5 = 1 corresponds to the spherical droplet, and we need to find < a < 1. 

In Fig. [2] the value of a that minimizes the energy, oq, is plotted as a function of the 
magnetic field in reduced units. The surface tension of a Meissner liquid contains a contribu- 
tion (Jsc coming from the interface between the superconductor and the normal state; this 
should be proportional to the difference between the coherence length and the penetration 
depth. There is also a contribution ai from the liquid-vapour interface. For the Meissner 
liquid reported in Ref. jsl, the surface tension is estimated to be of the order of a = 10~^ 
J/m^. For a 1 mm droplet, the applied magnetic flux corresponding to F = 1 is then of the 
order millitesla. Thus, fields reasonably below the critical field will induce a non-negligible 
deformation of the droplet. The smaller the droplet, the more it resists shape deformations. 

A Taylor expansion of the energy around the minimum allows to find the oscillation 
frequency of the quadrupole mode of the droplet. We equate the second order term in the 
expansion with a harmonic oscillation around the optimal shape. 
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FIG. 2: (color online) The droplet shape, characterized by the semiminor axis expressed in units 
ro of the spherical droplet with the same volume, is shown as a function of the magnetic field, 
expressed through F = ^H^rQ/(2a). In the inset, the oscillation frequency of the quadrupole mode 
(in units of loq = \l'ia j (2prQ)) is shown as a function of F. The limiting value for a spherical 
droplet is given. 

Here p is the density of the Meissner liquid. The dependence of the frequency on the magnetic 
field in reduced units, F, is shown in the inset of Fig. [2l The frequency is expressed in units 
uoq = ^J'iaj {2prQ). In the limit of spherical droplets {Ha — > 0) the frequency converges to 
uj/ujo = 4:^/2/5. For a ro =1 mm droplet with a = 10 ^ J/rn^ and p = 10^ kg/m^, the unit 
ujo corresponds to a frequency of 39 Hz. Increasing the magnetic field deforms the bubble, 
and also stiffens the oscillation frequency. 

III. SURFACE WAVES ON A MEISSNER FLUID 

The quadrupole oscillation mode of the droplet is one particular realization of surface 
waves. In this section, we investigate how the surface waves on an infinitely deep Meissner 
fluid are influenced by an applied magnetic field parallel to the surface and parallel to the 
propagation direction of the wave (a field normal to the propagation direction of the wave 
would be parallel to the wave fronts and would not induce magnetic charges on the surface - 
thus it would not alter the wave dispersion). Here, we consider small-amplitude waves. The 
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FIG. 3: (color online) The magnetization of a Meissner fluid with a surface wave is shown in 
panel (a); the applied field is parallel to the surface and to the wave propagation direction. The 
magnetization induces surface charges on the rising and descending wave slopes, as indicated. The 
shading (from dark blue - low to light green - high) indicates the concentration of magnetic charge. 
In panel (b) and (c), the demagnetization field and the total magnetic flux are shown, respectively. 

surface of the liquid is characterized by ({x) = a;sin(A;x) where the flat surface corresponds 
to the xy-plane {at z = 0), a is the wave amplitude, and A; = 27r/A is the wave number. The 
angle between the tangent to the surface and the horizontal plane is 

9{x) = arctan[C'(2;)] = ak cos{kx) + O(a^). (16) 
A. Demagnetization field 

When the surface of the fluid is flat and the applied field is parallel to the surface, 
the magnetization exactly cancels = HaG^ to achieve the Meissner state. However, 
when a wave is present, the magnetization will induce surface magnetic charge on the rising 
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and descending slopes of the wave, as illustrated in panel (a) of Fig. ([3]). The shading of 
the liquid in this figure also illustrates the concentration of this magnetic charge V ■ M, 
strongest at the slopes of the wave. As in a regular magnetized object, this gives rize to a 
demagnetization field H^, calculated below and shown in panel (b) of Fig. The resulting 
total flux is shown in panel (c) of that Fig. ([2]), and, as calculated below, follows the contour 
of the surface wave. 

To find the magnetic field H(a;, z) in this case, we decompose this field in a component H/f 
tangential to the surface and a component normal to the surface. For a small-amplitude 
wave, H± will be of order a and H// ^ \Ha\ up to second order in a. The demagnetization 
field Hrf = H — Hae^ at the surface of the wave can then be written as 

= {H,/ cos [e{x)] - H^{x) sin [e{x)] - H,] 

+ {H/i?>m[e{x)]+ H^{x)cos[e{x)]]e,, (17) 
or, expanding with respect to a, 

H, [x,C(x)] 



-Ha-a'^k'^ cos^{kx) — H^{x)ak cos{kx) 



(18) 



+ [Haakcos{kx) + H±{x)] + 0{a^). 

For the calculation of the energy associated with the magnetization, we need to find the 
demagnetization field inside the fluid. As a trial solution, we use 

iif{x, z) = h(a;)e'=(^-'^(^)). (19) 

That is, we assume an exponential decay into the fluid with the wave length of the surface 
wave as a characteristic scale. In the direction of propagation of the wave, we need to find 
h{x) = hx{x)ex + hz{x)ez. The ansatz f|T9|) allows to satisfy V x H™ = 0, which reduces to 



dhzjx) 
dx 



— hz{x)ak^ cos{kx) = kh^^x), 



(20) 



with boundary conditions h{x) = H(i[x,C(a;)] and h(x) = h{x + A). This can be solved 
straightforwardly, yielding 

Hrf [x, z < C{x)] = Ha^— cos^(/cx) 



X exp{/c [z — a sm{kx)]}ex 



(21) 
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for the demagnetization field inside tlie liquid. 



B. Magnetic energy 

Consider a surface of fluid of length Ly in the y-direction and = p\ {p E N) in 
the x-direction. The energy difference Emag between the unmagnetized (Ha = 0) and the 
magnetized case can be written as 

E^^g = ^ /" Ha ■ H,dh + ^ / nldh. (22) 



2 J 2 _ 

z<((x) 

The second term represents the energy of swichting on the applied magnetic field, and does 
not depend on the presence of a surface wave. The first term evaluates to 



Emag = ^ J dx COS (fcx) 



2 



k LrL., a 



2 



4 2 

Thus, the magnetic energy per unit surface required to establish the Meissner state when 
there is a surface wave with wave number k and amplitude a is (to order a^) given by 

Em = ^ka\ (23) 

C. Wave dispersion 

To find the dispersion relation of a surface wave on an infinitely deep Meissner fluid, we 
follow the Hamiltonian procedure outlined in Ref. [lOi]. The kinetic energy per unit surface, 
associated with the surface wave is 

Ek = \^a'. (24) 

When the flat surface is deformed, restoring forces tend to pull it flat again. These restoring 
forces can be related to the surface tension energy, per unit surface: 



= ^ak^a'^, (25) 
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and to the gravitational energy, per unit surface: 

Eg = \p9a'. (26) 

In a Meissner liquid, subjected to a magnetic field parallel to the surface wave propagation 
direction, also the magnetic energy ( l23l) gives rise to a restoring force. The (classical) 
Hamiltonian associated with the surface wave of amplitude a is then given by the sum of 
((21,(I251),(I26D and ((231) : 

This expression is valid for small-amphtude oscillations, and leads to a dispersion relation 

^ik) = J + ^-+gk. 28 

VP 2 p 

The effect of the magnetic field dominates when jj^Ha » ^/Jipgjk and ^Ha » \^JIak. 
Taking again a = 10~^ J/m^ and p = 10^ kg/m^, and a magnetic field of fiHa = 100 G, we 
find that the Meissner contribution to the dispersion dominates for 10^ m~^ ^ k < 10^ m~^. 
The experimental setup of Ref. jj] would allow to probe this range of wave lengths. Surface 
waves can be detected by reflecting a laser beam off the surface of the Meissner fluid and 
noting the displacement of the laser spot as a function of time and space. 



IV. CONCLUSION 



Fluids with remarkable magnetic response properties, such as ferrofluids, have sparked 
a lot of interest. Here, we investigate how a fluid superconductor, a Meissner fluid, would 
react to an applied magnetic field. Even though such Meissner f uids are not yet accessible, 
candidate realizations from theory [2!] and recent experiments js] can be found. 

When such a fluid superconducting material is placed in an applied magnetic field, it will 
flow to adapt its shape and minimize the energy required to expell the magnetic flux. In this 
contribution, we focused on a single droplet of Meissner fluid and derived the deformation 
of such a droplet due to an applied magnetic field as well as the associated quadrupole 
oscillation frequency. Inversely, imposing a hydrodynamic flow, for example by generating 
a surface wave, will alter the energetics and thus the dispersion of the wave. We calculated 
how the dispersion of the surface wave is changed due to the presence of a magnetic field, 
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applied parallel to the surface and to the propagation direction of the wave. The change in 
the dispersion is predicted to be relevant for the wave length range thought to be achievable 
in a setup such as ^ . 
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